Thermodynamics of a heavy ion-irradiated superconductor: the zero-field transition 
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Specific heat measurements show that the introduction of amorphous columnar defects consid- 
erably affects the transition from the normal to the superconducting state in zero magnetic field. 
Experimental results are compared to numerical simulations of the 3D XY model for both the pure 
system and the system containing random columnar disorder. The numerics reproduce the salient 
features of experiment, showing in particular that the specific heat peak changes from cusp-like 
to smoothly rounded when columnar defects are added. By considering the specific heat critical 
exponent a, we argue that such behavior is consistent with recent numerical work^ showing that 
the introduction of columnar defects changes the universality class of the transition. 



PACS numbers: 74.25.Bt,74.40.+k,68.35.Rh 

An enormous amount of attention has been paid in 
recent years to the effect of amorphous columnar de- 
fects on the superconducting transition in a magnetic 
K(T]i— — — — Much less work has been done on the 
transition in zero magnetic field, the only report to our 
knowledge being the measurement of the microwave con- 
ductance transverse to the columnar defects^ Neverthe- 
less, the transition in zero magnetic field merits attention 
in its own right £ It is expected that disorder will be a rel- 
evant perturbation, and change the universality class of a 
phase transition, whenever 2 — d*v > (modified Harris 
criterion^) . Here d* is the number of dimensions in which 
the system is disordered, and v is the usual correlation 
length critical exponent. In the absence of columnar de- 
fects, we expect the superconducting phase transition to 
fall in the universality class of the three dimensional (3D) 
XY model, so that v = 0.6717(1) J£ In the case of random 
point disorder, d* = 3, so that 2 - d*v « -0.015 < 0, 
and the disorder is irrelevant. In the case of columnar 
disorder, however, d* = 2, so that 2 - d*v w 0.66 > 0, 
and disorder should be relevant and drive the system to 
a new universality class. Note that the stability of the 
new disordered critical point with respect to the modified 
Harris criterion requires that the new correlation length 
critical exponent satisfies v > l£ Recent simulations of a 



columnar disordered XY model in Ref. d supported such 
expectations, finding a phase transition with anisotropic 
scaling and a value for the critical exponent v w 1.2. 

It is the purpose of this report to test these ideas on 
a real superconducting system. For this, we have chosen 
to measure the specific heat of optimally doped, single 
crystalline YBa2Cu30y_5, both without and with colum- 
nar defects. It has been shown previousl y 11 ^ 2 that the 
specific heat at the superconducting transition of pris- 
tine YBa2Cri307_5 is consistent with that of the 3D 
XY model, with a small negative specific heat exponent 
a w — 0.015J£ It turns out that this exponent is con- 
siderably modified by the introduction of the columnar 
defects. 

Experiments were done on two YBa2Cu30y_5 single 
crystals, cut from the same piece. The original crystal 
was grown by the flux method in Au crucibles, and sub- 
sequently annealed in oxygen in Pt tubesJ^ One crystal, 
which contained a single family of twin boundaries sep- 
arated by a distance of approximately 10 fj,m, was irra- 
diated with 5.8 GeV Pb ions at the Grand Accelerateur 
National d'lons Lourds (GANIL) in Caen, France, to a 
fluence of 1 x 10 11 ions cm" 2 . The ion beam was directed 
parallel to the c-axis; each ion impact created an amor- 
phous columnar track of radius ~ 3.5 nm. The second, 
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untwinncd crystal, was not irradiated, but kept as the 
pristine reference sample. Specific heat measurements 
have been performed in the absence of an applied mag- 
netic field, using the same measurement technique em- 
ployed in Ref . [lj. 

In Fig. QJi we show our raw data for the specific heat of 
both the pristine and irradiated crystals. The tempera- 
ture axis has been scaled by the value of the temperature 
Tpcak at the specific heat peak, in order to better com- 
pare the two samples (T peak = 93.1 K for the pristine 
sample, while T pGa k = 92.1 K for the irradiated sam- 
ple). We see clear specific heat anomalies, signatures of 
the superconducting phase transition, superimposed on 
a smooth increasing background. In the pristine sam- 
ple, the amplitude of the specific heat anomaly was of 
the order of 4% of the total specific heat, attesting to 
its very high quality. Fitting the smooth background 
to a cubic polynomial (dashed line in Fig. [TJl) we sub- 
tract this background from the data in Fig. [TJb, in order 
to better emphasize the shape of the anomaly. For the 
pristine sample, we see the typical "A" -cusp shape ex- 
pected for the H = superconducting transition in the 
presence of strong thermal fluctuations. The introduc- 
tion of the amorphous columnar defects reduces the ab- 
solute temperature at which the specific heat maximum 
occurs. Note that a lowering of the critical temperature 
after heavy ion irradiation may occur as a result of "self- 
doping" of the intercolumn space by O ions expelled from 
the tracks j 15 ' 16 however, no such effect was reported for 
YBa2Cu307_,5. Another possibility is that the columns 
reduce the average T c at which long range superconduct- 
ing order can set ini 17 i 18 The resulting specific heat-curve 
after heavy-ion bombardment shows notable differences 
with respect to the curve before irradiation. Most specif- 
ically, we find that the shape of the maximum is now 
smoothly rounded rather than the sharp cusp seen in the 
pristine sample. Following the suggestion^ that the in- 
troduction of columnar defects changes the universality 
class of the superconducting transition in zero field, we 
propose that the associated change of critical exponents 
is at the origin of the markedly different shape of the 
specific heat peak before and after irradiation. 

The presence of columnar defects implies that, even in 
zero magnetic field, critical scaling of physical quantities 
may be anisotropic: the correlation length parallel to the 
columns, £ z , diverges as a different power of the reduced 
temperature t = (T — T c )/T c than does the correlation 
length in the transverse direction, This defines the 
anisotropy exponent £, 
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Defining the correlation length exponent v in the usual 
way, £_l ~ \t\~ v , the singular part of the free energy 
density will scale as 



f( T ) ~ ~ £± 2 ~ c 
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FIG. 1: (a) Raw specific heat data (in arbitrary units) taken 
on a pristine YBa2Cu307_<s crystal (•), and on a crystal with 
a columnar defect density rid = 1 x 10 11 cm -2 , i.e. a matching 
field B^, = $o^d = 2 T (x). In order to compare the two 
data sets, the temperature has been rescaled by the value 
Tpeak at which the specific heat peak occurs. The dashed line 
represents the smooth background contribution to the specific 
heat, which is estimated using a third order polynomial fit. 
This background is subtracted in (b) in order to emphasize 
the shape of the specific heat peak due to the superconducting 
transition. 



will scale as 
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As a consequence, the specific heat per unit volume, c, 



Thus, the specific heat exponent in the anisotropic case 
is a = 2 — v{2 + Q. For the pure system, i.e. a su- 
perconductor without columnar defects, the anisotropy 
exponent £ = 1, and v > |. The specific heat exponent 
a is small and slightly negative. In contrast, the calcula- 
tion for the XY model of a superconductor with columnar 
defects of Ref. HI found values C ~ 1.3 and v « 1.2. In this 
case, a ~ —2.0, and is thus much more strongly negative. 
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FIG. 2: Monte-Carlo calculation of the specific heat using 
the 3D XY model, for a pure (pristine) superconductor (o), 
and a superconductor containing columnar disorder. For the 
pure case, a system of size 40 3 was used, while for the dis- 
ordered case the calculation was done for two system sizes, 
with L = 40 (x) and L = 20 (□) grid points in the direction 
perpendicular to the anisotropy axis, respectively; parallel to 
the anisotropy axis, L z — 0.5L 1 ' 3 . The temperature axis has 
been rescaled by the value T pea k at which the specific heat 
peak occurs. 

Because a is negative for both cases, the specific heat 
does not diverge at the transition. However, an inter- 
esting difference is seen if we consider the temperature 
derivative of the specific heat, 



For the pure superconductor, — a - 1 « —0.985 is nega- 
tive. Therefore, the slope of the specific heat diverges at 
T c , giving rise to the familiar cusp observed in Fig. [1] for 
the pristine sample. For the superconductor with colum- 
nar defects however, —a — 1 w 1.0 is positive, so the slope 
of c does not diverge. There is no sharp cusp, as is indeed 
observed experimentally^ 

The above discussion considers only the singular part 
of the specific heat that comes from the large length 
scale critical fluctuations. It disregards the smooth non- 
singular contribution that comes from the non-critical 
short length scale fluctuations. Most likely, this non- 
singular part of the specific heat has some non-zero tem- 
perature derivative at T c . For the unirradiated supercon- 
ductor, this temperature derivative is much smaller than 
the diverging slope of the singular part, and can there- 
fore be disregarded. Thus, the non-singular contribution 
is unlikely to affect the shape of the c curve very much. 
In the presence of columnar defects, however, the slope 



of the singular part vanishes at T c . Then, the temper- 
ature derivative of the specific heat at T c is determined 
by the non-singular part. Since the slope of c is non-zero 
and smooth at T c , the maximum of c is no longer located 
at the critical temperature, and the shape of the peak is 
now a smoothly curved maximum. As we discuss below, 
this feature too is found in the experimental data. 

To illustrate the comparison further, we have carried 
out Monte Carlo simulations to numerically compute c 
for the pure system, and also for the system with colum- 
nar disorder, using the same XY model and random dis- 
tribution for the columnar disorder as defined in Sec. IIA 
of Ref. H For the pure case, we use a system size of 40 3 
grid sites. For the case with columnar disorder we used 
two different system sizes. Taking L = L x = L y as the 
system length transverse to the direction of the columnar 
defects, and using L z = 0.5L^ (£ = 1.3) consistent with 
anisotropic scaling, we considered sizes L = 20, averaged 
over 125 different realizations of the random disorder, 
and L = 40, averaged over 66 different realizations of 
disorder. Our results are shown in Fig. O It is clearly 
seen that the pure case displays a peak with the famil- 
iar "A" -cusp, whereas the case with columnar disorder 
displays a smoothly rounded peak. We observe essen- 
tially no size dependence in our results for the case with 
columnar disorder. Our numerical results thus repro- 
duce the main qualitative features found experimentally 
for single-crystalline Yl^CuaOy-^. Difference in detail 
between the experimental and model curves for the case 
with columnar defects most likely results from the sim- 
plicity of the numerical model as well as differences in the 
effective strength of the disorder. For the simulations, a 
particularly strong disorder was chosen so as to reach the 
asymptotic limit even with the small system sizes under 
scrutiny. Presumably, this explains why the model curve 
is even rounder than the experimental one. 

Another feature of the numerical results is that, in 
agreement with the argument above, the system with 
columnar disorder is found to have a specific heat peak 
that lies below the critical temperature, with (T c — 
Tpeak) /T c = 0.07. In order to compare this result with ex- 
periment, the critical temperature of the YB^CuaOy-a 
crystals is deduced from the functional dependence of 
the thermodynamic properties of the superconductor on 
a relevant "scaling" parameter. In Ref. [20j , it was shown 
that, for magnetic fields B larger than 1 T, the field 
and temperature dependence of the magnetization of 
the crystals under study can be described by a unique 
(Ginzburg-Landau Lowest Landau Level, GL-LLL) re- 
lation, M/($b) 2 / 3 cx F(Q), with Q = (1 - b)(l - 
z? 2 ) 1/3 (^)~ 2/3 , where = T/T c and b = B/B c2 (T); here 
T c is the zero field critical temperature and B C 2(T) the 
upper critical field. For the unirradiated crystal, this re- 
lation gives a critical temperature T c = 93.1 K, which 
coincides to within 0.1 K with the position of the zero 
field specific heat maximum. For the irradiated crys- 
tal, a similar analysis gives, within the accuracy of the 
fit, the same critical temperature T c = 93.1 K, which 
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now lies above the observed specific heat maximum at 
92.1 K. The experimental data thus show a separation 
(T c — T pea k)/T c = 0.01, in the same direction as the 
numerical results, but smaller in magnitude, again pre- 
sumably due to the larger disorder strength used in the 
numerical model calculations. 

In summary, we have measured the specific heat 
of pristine and heavy-ion irradiated single crystalline 
YF^CujOt-^ in zero applied magnetic field. The re- 
sults were compared to Monte Carlo simulations of the 
3D XY model for both the pure case and the case with 
columnar disorder. Both experiment and numerics show 
a drastic influence of the columnar defects on the shape 
of the specific heat anomaly at the superconducting-to- 



normal transition. The overall features of the specific 
heat anomaly are well explained by the critical expo- 
nents obtained from the numerical calculations, suggest- 
ing that the introduction of the columnar defects does 
indeed change the universality class of the transition. 
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